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Abstract 


The  atomic  structure  and  energetics  of  the  clean  (100)-(2xl),  (100)-(2xl):H  roonohydride, 
(l(X))-(lxl):2H  full  dihydride,  (100)-(3xl):1.33H  intermediate  dihydride,  and  other  intermediate 
hydride  surfaces  obtained  by  local  removal  of  H  atoms  from  or  addition  of  H  atoms  to  the  (2x1) 


monohydiide,  are  investigated  by  molecular  mechanics  (MM3).  The  monohydride  phase  is  found 
to  be  the  most  stable  thermodynamically  and  is  predicted  to  be  the  dominant  phase  under  chemical 
vapor  deposition  conditions.  Gas-surface  reactions  relevant  to  diamond  chemical  vapor  deposition 
involving  hydrogen  and  the  diamond  (l(X))-(2xl):H  surface  are  discussed  in  light  of  the  derived 
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I.  INTRODUCTION 

In  the  past  few  years,  our  understanding  of  the  low-pressure  growth  of  diamond  by  chemical 
vapor  deposition  (CVD)  has  been  greatly  enhanced  due  to  intensive  research  efforts.  ^  Much  has 
been  learned  about  the  gas  phase  chemistry  during  diamond  CVD^  and  the  identity  of  the  primary 
growth  precursor(s),3  particularly  for  hot-filament  CVD.  Correspondingly,  many  mechanistic 
models  based  on  the  observed  gas-phase  species,  such  as  methyl  radical'^  and  acetylene,^  have  been 
proposed  to  account  for  diamond  growth  at  the  molecular  level.  In  contrast,  very  little  is  known 
about  the  chemistry  on  the  diamond  surfaces,  especially  under  growth  conditions.  In  most  of  the 
mechanistic  models  thus  far  proposed,  very  simple  assumptions  about  the  atomic  structure  of  the 
surfaces  under  grov/th  conditions  are  made  and  the  kinetics  and  thermodynamics  of  surface 
reactions  are  estimated  fiom  gas-phase  reactions.  It  is  clear  that  considerably  more  experimental 
and  theoretical  work  will  be  required  to  understand  diamond  growth  at  the  molecular  level. 

Among  the  three  low-index  surfaces  of  diamond,  i.e.,  (1 1 1),  (1 10),  and  (100),  the  diamond 
(100)  orientation  is  uttique  because  nearly-atomically-smooth  surfaces,  critical  to  comparison 
between  experiment  and  theory,  have  been  obtained  only  on  the  (100)  face  by  homoepitaxial  CVD 
growth.  3ti.6-8  Moreover,  each  of  carbon  atoms  on  the  ideal  (100)  surface  has  two  dangling  bonds 
(unpaired  electrons)  while  on  ideal  (110)  and  (1 1 1)  surfaces  each  surface  carbon  atom  has  only  one 
dangling  bond,  thus  rendering  rich  and  yet,  up  to  now,  very  confusing  chemistry  which  provides  a 
challenge  for  both  experiment  and  theoiy.  Experin^ntal  studies  of  diamond  (100)  have  been  carried 
out  by  low-energy  electron  diffraction  (LEED),^'!!  temperature-programmed  desorption 
(TPD),10*ll  photon-stimulated  ion  desorption  (PSID),12  electron-stimulated  desorption  time  of 
flight  (ESD-TOF)  spectroscopy,  photoemission,  and  Rutherford  backscattering  spectroscopy 
(RBS).^3  But  many  important  questions,  such  as  the  surface  hydrogen  coverage,  assignment  of 
TPD  desorption  peaks,  and  stability  of  hydride  phases,  are  still  issues  of  debate.  On  the  theoretical 
side,  various  authors  have  utilized  semiempiricall^l^  and  empirical  approaches  1^*19  to  investigate 
the  properties  of  clean  and  hydrogenated  diamond  (100).  But  substantial  discrepancies  remain,  in 
both  the  dimer  bond  length  (±  0.1  A)  and  in  the  energetics,  and  a  satisfactory  description  of  the 
chemisorption  and  reaction  of  hydrogen  on  diamond  (100)  is  not  yet  possible. 

We  recently  demonstrated  the  usefulness  of  the  empirical  molecular  mechanics  (MM3)  method, 
developed  by  AUinger  and  co-workers  over  the  past  two  decades,2021  in  studying  the  surface 
chemistry  of  diamond  (lOO).^^  closely-related  approach  has  been  taken  by  Harris  and  Belton  in 

calculations  on  diamond  (1 1 1)  and  (1 10)  surfaces.^^  Molecular  mechanics  is  arguably  the  most 
successful  and  certainly  the  most  widely  used  of  the  empirical  potentials  available,  and  several 
versions  are  commercially  available.  In  the  present  paper,  further  thermochemical  analysis  is 
performed  on  the  previous  results,  and  several  additional  hydride  stmctures  are  considered.  Recent 
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near-atomic-resolution  images  of  CVD-grown  diamond  (100),  obtained  by  scanning  tunneling 
microscopy^  and  atomic  force  microscopy,®  show  evidence  of  a  (2x1)  reconstruction,  presumably 
associated  with  the  so-called  monohydride  structure  with  one  hydrogen  atom  per  surface  carbon 
atom  and  a  dimer  bond  between  pairs  of  surface  carbon  atoms7'l0>^'^'19«22  Because  of  the  likely 
importance  of  the  monohydride  surface  in  CVD  growth,  H  atom  abstraction  and  insertion  reactions 
on  this  particular  surface  are  studied  by  molecular  mechanics  as  well.  Several  gas-surface  reactions 
important  in  CVD  growth  are  discussed  in  light  of  the  derived  energetics. 

II.  THEORETICAL  METHOD  AND  RESULTS 

The  calculation  method  used  to  calculate  the  structures  and  enthalpies  of  formation  of  the 
diamond  (100)-(2xl),  (l()0)-(2xl):H,  (l()0)-(3xl):1.33H,  and  (100)-(lxl):2H  surfaces  has  been 
described  in  detail  previously.^  Die  third-generation  MM3  force  field,  with  parameters  for 
saturated,  unsaturated,  and  conjugated  hydrocarbons,  has  a  demonstrated  high  degree  of  accuracy 
(bond  lengths  «  ±0.01  A,  bond  angles  between  atoms  other  than  hydrogen  »  ±1%  torsional  angles 
«  ±4*,  heats  of  formation  »  ±1  kcal/mol)  for  small,  large,  and  highly  strained  molecules  and  bulk 
diamond  as  well.20.21  mM3  should  be  applicable  to  the  description  of  saturated,  unsaturated,  and 
conjugated  hydrocarbon  species  on  any  crystal  face  of  diamond  as  long  as  the  bond  lengths,  bond 
angles,  and  distances  between  nonbonded  atoms  are  within  the  range  of  values  in  structures  for 
which  MM3  has  demonstrated  accuracy.  MM3  paranxters  for  radicals  are  tentative, making 
calculations  with  open-shell  species  more  uncertain,  and  MM3  cannot  describe  surface  species  with 
bonding  configurations  that  have  not  been  parameterized  in  molecules. 

Briefly,  atomic  positions  are  determined  by  minimizing  an  empirical  potential  energy  function, 
the  so-called  steric  energy  E,202l  subject  to  periodic  boundary  conditions'^  and  a  substrate  lattice 
constant  fixed  at  the  value  of  bulk  diamond.  E  is  given  by 

E  =  'ZEs  +  XEe  +  XEo)  +  XEs6  +  XEo3s  +  +  Z^vdW  +  Z^dpl  (1) 

where  Es  is  the  bond  stretching  energy,  Eg  is  the  angle  bending  energy,  Eais  the  torsional  energy; 
EsS  is  the  stretch-bend  energy,  Eos  is  the  torsion-stretch  energy,  Eg^  is  the  bend-bend  coupling 
energy;  EvdW  is  the  van  der  Waals  repulsion  energy  between  nonbonded  atoms,  and  Edpl  is  the 
dipole-dipole  interaction  energy  .20^1  The  enthalpy  of  formation  is  calculated  by  adding  bond 
enthalpies  (Af/bond)  and  functional-group  correction  terms  (AZ/gtruct)  to  the  minimized  steric  energy: 

AHf  =  A//bond  +  A//smici  +  E.  (2) 

For  completeness  and  comparison  to  the  new  results,  scale  drawings  of  the  structures  of  the 
clean  (100)-(2xl),  (100)-(2xl);H  monohydride,  (100)-(lxl):2H  full  dihydride,  and  (100)- 
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(3xl):1.33H  partial  dihydride  surfaces  are  presented  in  Figure  1,  with  selected  bond  lengths  (A) 
and  angles  given  explicitly.  The  calculated  energetics  of  these  structures  are  described  by  the 
following  enthalpies  of  reaction  at  298  K,  expressed  with  respect  to  the  (2xl);H  monohydride:^^ 

C(100)-(2xl):H  C(100)-(2xl)  +  H2  AH 29%  =  +46.7  kcal/mol  (3) 

C(100)-(2xl):H  +  H2  ^  2  C(100)-(lxl):2H  A// 298  = +49.2  kcal/mol  (4) 

3  C(100H2xl):H  +  H2  2  C(100)-(3xl):1.33H  AH 29%=  -15.6  kcal/mol  (5) 

We  have  also  calculated  temperature-dependent  enthalpies,  entropies,  and  fiee  energies  of 
reaction  for  Eqs.  (3-5),  enabling  estimation  of  equilibrium  concentrations  under  CVD  conditions: 


AH'(J)  =  +  1  AC'p  dr 

(6) 

^298 

AS'(J)  =  AS  ^98  +  I  AC'pdlnr 

(7) 

J79S 

AG’  =  AH'  -  TAS' 

(8) 

Typically,  entropy  and  heat  capacity  changes  are  only  weakly  dependent  on  the  steric  energy,  and 
we  have  estimated  AS  298  gas-surface  reactions  (3-5)  at  various  ten^ratures  using 

the  group  additivity  method.^^  The  results  are  summarized  in  Table  I.  AH'  is  only  weakly 
temperature-dependent,  and  therefore  small  errors  in  AC'p  are  inconsequential.  The  largest  contri¬ 
bution  to  AS  298  in  Eqs.  (3-5)  comes  from  the  lost  translational  entropy  of  gas-phase  H2,  and 
therefore  the  group  additivity  estimates  should  be  fairly  accurate.  In  cases  whwe  the  steric  energy  is 
not  small,  e.g.,  the  (lxl);2H  dihydride,  the  group  additivity  estimates  are  much  less  accurate.  The 
entropy  associated  with  opening  polycyclic  ring  systems  is  difficult  to  estimate  using  group 
additivity,  as  pointed  out  by  Harris  et  and  the  AS  estimates  for  reactions  (4)  and  (5)  may  be 
uncertain  by  5-10  cal  mol'^  K'l. 

We  believe,  based  on  the  currently  available  experimental  evidence^’^'^^  and  the  present 
theoretical  work,  that  the  (2xl):H  monohydride  surface  (Fig.  1(b))  is  the  predominant  species 
under  typical  hot-filament  or  microwave  plasma  CVD  growth  conditions.  Therefore,  we  have 
performed  MM3  calculations  on  surfaces  which  are  related  to  the  monohydride  by  removal  or 
addition  of  hydrogen  atoms  due  to  abstraction  by  or  insertion  of  gas-phase  H  atoms,  respectively. 

In  contrast  to  our  previous  study ,22  we  have  chosen  the  (l(X))-(2xl):H  monohydride  surface  as  the 
reference  state,  so  that  the  surface  steric  energy  is  defined  as  the  difference  between  the  steric 
energy  of  species  obtained  by  removal  or  addition  of  H  atoms  on  the  monohydride  and  the  stoic 
energy  of  the  monohydride  itself. 


The  slab  chosen  for  the  initial  and  modified  diamond  (l(X))-(2xl):H  monohydride  surfaces  was 
5  layers  thick,  consisting  of  3x4  (2x1)  unit  cells  with  a  total  of  123  carbon  atoms.^^  To  simplify 
the  minimization  of  the  surface  steiic  energy,  we  allowed  only  the  surface  atoms  (C  and  H)  in  the 
unit  cell  which  either  loses  or  gains  H  atoms  to  relax  and  fixed  the  remainder  of  the  atoms  in  their 
monohydride  surface  positions.  The  surface  energy  calculated  in  this  way  is  estimated  to  be 
approximately  2-3  kcal/mol  less  accurate  than  that  obtained  by  complete  relaxation  of  the  surface 
atoms.  An  additional  uncertainty  of  about  ±3  kcal/mol  for  radical  surface  species  is  estimated. 

The  surface  species  formed  by  simple  removals  or  additions  of  hydrogen  atoms  on  the  (100)- 
(2xl):H  monohydride  surface  are  summarized  in  Eqs.  (9-12)  below,  together  with  a  shorthand 
notation  (“R”  denotes  a  backbond  to  second-layer  lattice  atoms)  and  the  calculated  enthalpies  of 
formation  from  the  monohydride. 


H  H 

(SR2CH-CHR2) 

R2CH-CHR2 

R2CH-CHR2  +  jH2 

R2CH-CHR2  +  H2  -> 


H 

(sR2CH-CR2)  +  7H2 
(=  R2C=CR2)  -1-  H2 

H  H  H 

(SR2CH2CHR2) 

H  HH  H 

Y  Y  (=R2CH2CH2R2) 


A// 298  =  39.3  kcal/mol  (9) 
AN  298  =  44.6  kcal/mol  (10) 

AH  ^98  =  37.6  kcal/mol  (11) 

AH  ^98  =  12.6  kcal/mol  (12) 


In  reaction  (9),  a  single  hydrogen  atom  is  removed  from  a  monohydride  surface  dimer,  creating  a 
radical  site,  while  both  hydrogen  atoms  are  removed  in  reaction  (10),  forming  a  strained  dimer  it 
bond.  In  reaction  (1 1)  the  dimer  bond  is  broken  and  a  hydrogen  atom  is  added  to  one  of  the  carbon 
atoms,  while  reaction  (12)  represents  a  hydrogenolysis  of  the  dimer  C-C  bond.  In  each  of  these 
stmctures,  all  surrounding  dimers  remain  in  the  (2x1  ):H  monohydride  configuration.  ITie 
calculated  atomic  structures  of  R2CH-  C  R2,  R2C=CR2,  R2CH2  C  HR2,  and  R2CH2CH2R2  are 
given  in  Fig.  2(a-d),  respectively,  and  the  displacements  of  atoms  from  their  ideal  bulk-taminated 
positions  are  presented  in  Table  II. 


III.  DISCUSSION 

A  detailed  discussion  of  the  structure  and  energetics  of  the  clean  (l()0)-(2xl),  (l()0)-(2xl):H 
monohydride,  (100)-(lxl):2H  full  dihydride,  and  (l()0)-(3xl):1.33H  partial  dihydride  surfaces 
(Fig.  1)  has  been  presented  elsewhere.^^  Briefly,  the  structures  are  in  reasonable  agreement  overall 
(±  0.1  A  in  bond  lengths)  with  the  recent  semiempirical  results  of  Verwoerd^^  and  of  Mehandru  and 


Andersonl^  and  the  empirical  results  of  Brenner,!^  with  larger  discrepancies  in  the  energetics.  The 
demonstrated  performance  of  MM3  for  a  wide  range  of  hydrocarbons^O  suggests  that  our  results  are 
the  most  accurate  of  the  existing  surface  calculations  for  structures  with  bond  lengths  and  bond 
angles  in  the  range  known  to  be  well  described  by  MM3  (a  similar  limitation  applies  to 
semiempirical  theories  and  other  empirical  potentials).  The  (lxl):2H  full  dihydride  and,  to  a  lesser 
extent,  the  clean  (100)-(2xl)  structures  have  parameters  (nonbonded  H-H  distance  and  C=C  bond 
length,  respectively)  that  lie  outside  the  established  range,  and  our  calculations  are  less  reliable  than 
for  the  (2xl):H  and  (3xl):1.33H  structures.  The  very  recent  semiempirical  calculations  of  Zheng 
and  Smith^^  show  substantial  discrepancies  from  the  aforementioned  studies,l^*^6-*9,22  including 
significantly  shorter  C-C  dimer  bond  lengths  and  a  prediction  of  a  highly  stable  (100)-(lxl):2H 
dihydride  (A//*  at  0  K  =  -86.6  kcal/mol  for  reaction  (4)^^). 

It  is  apparent  fiom  Table  I  that  among  the  three  possible  diamond  surface  hydride  phases,  the 
monohydride  phase  is  the  most  stable  over  the  temperature  range  of  298-1500  K,  followed  in 
stability  by  the  clean  surface,  which  suggests  that  the  monohydride  is  the  dominant  surface  species 
at  the  typical  growth  tempertatures  of  1073-1273  K.1  Recent  STM,^  AFM,^  andUEED^^ 
investigations  on  homoepitaxial  diamond  (100)  films  transferred  and/or  analyzed  in  air  have  clearly 
observed  a  (2x1)  reconstruction  with  dimer  rows  and  steps,  which  each  set  of  authors  has  attributed 
to  the  monohydride.7.8.11  While  observation  of  a  surface  structure  after  growth  and  air  exposure 
does  not  prove  that  the  structure  was  present  under  growth  conditions,  it  is  strongly  suggestive, 
particularly  when  that  structure  is  a  stable  one  and  its  appearance  is  insensitive  to  the  details  of 
growth  and  reactor  shutdown  conditions,  as  evidenced  by  its  observation  by  several  research 
groups.  Since  less-stable  structures  like  the  (3xl):1.33H  and  (lxl):2H  dihydrides  arc  favored  at 
lower  temperatures,  if  they  did  exist  under  growth  conditions  they  would  likely  be  stable  as  growth 
was  quenched  and  the  substrate  was  cooled.  The  fact  that  they  arc  not  observed  provides  indirect 
evidence  that  they  also  do  not  exist  in  large  concentrations  under  growth  conditions  and  that  the 
(2xl):H  monohydiide  predominates,  with  clean  (2x1)  dimers  perhaps  also  being  important. 

An  important  role  for  dihydride  phases  on  diamond  (100)  has  been  suggested  by  several 
authors.^^40,17  Dihydride  species  are  presumably  generated  in  a  local,  stepwise  manner  from  the 
monohydride  by  reaction  with  atomic  hydrogen,  so  it  is  appropriate  to  first  consider  the  local 
modifications  of  the  (2xl):H  monohydride  illustrated  in  Fig.  2.  The  energetics  of  these  intercon¬ 
versions  follow  directly  from  Eqs.  (9-12)  and  the  heat  of  formation  of  H  (52.10  kcal/mol^^): 


R2CH-CHR2 

+  H 

R2CH-CR2  +  H2 

A// 298  =  -12.8  kcal/mol 

(13) 

R2CH-CR2 

+  H 

R2C=CR2  +  H2 

A// 298  =  -46.8  kcal/mol 

(14) 

R2CH-CHR2 

+  H 

-> 

R2CH2CHR2 

A// 298  =  -14.5  kcal/mol 

(15) 

R2CH2CHR2  +  H  R2CH2CH2R2  A//298  =  -77.1  kcal/mol  (16) 

When  a  H  atom  on  the  (2x1)  monohydride  surface  is  abstracted  by  a  gas-phase  H  atom,  as  in 
reaction  (13),  a  radical  site  is  formed  on  the  surface.  The  radical  C  atom  moves  towards  the  dimer 
center  and  slightly  upward  from  the  (2x1  ):H  position  towards  a  more  favorable  sp2  geometry, 
while  the  neighboring  C  on  the  same  dimer  moves  away  from  the  dimer  center.  The  change  in 
enthalpy  associated  with  abstraction  of  H  from  the  monohydride  surface,  -12.8  kcal/mol,  is  similar 
to  but  somewhat  more  negative  than  the  value  associated  with  abstraction  of  the  tertiary  hydrogen 
atom  from  isobutane  to  form  a  rm-butyl  radical,  -8.4  kcal/mol.^  The  difference  is  likely  due  to 
release  of  some  of  the  strain  in  the  monohydride.  Brenner,  using  an  empirical  Tersoff-type 
potential,  obtained  - 10.0  kcal/mol  for  the  same  abstraction  reaction.  Since  his  potential  was 
constmcted  using  atomization  energies  of  hydrocarbon  molecules  at  0  K,  a  temperature  correction 
should  be  made  to  his  value  in  order  to  make  a  direct  comparison  with  our  value.  Brenner  also 
found  that  the  enthalpy  of  abstraction  of  one  H  atom  from  the  diamond  (1 1  l)-(lxl):H  monohydride 
surface  was  nearly  the  same  (A//  =  -11  kcal/mol)  as  on  the  (100)  face.  In  both  cases,  the  H  atom  is 
attached  to  a  tertiary  carbon  atom  and  therefore  the  abstraction  enthalpies  should  be  similar.  Both  the 
present  results  and  those  of  Brennerl^  support  the  applicability  of  gas  phase  kinetic  and 
thermodynamic  data  to  simulation  of  diamond  growth  kinetics^l’»'W*73.27  also  indicate  that  the 

energetics  can  easily  differ  by  several  kcal/mol  between  analogous  molecular  and  surface  species. 

Both  the  structure  and  energetics  of  the  R2C=CR2  species  formed  by  removal  of  two  hydrogen 
atoms  from  the  (2xl);H  monohydride  are  quite  similar  to  those  of  the  clean  surface,  indicating  that 
our  present  assumption  of  no  subsurface  relaxation  is  not  too  severe.  For  comparison,  we  calculate 
a  dimer  bond  length  of  1.48  A  in  R2C=CR2  (surrounded  by  R2CH-CHR2  in  adjacent  unit  cells) 
versus  1.46  A  on  the  clean  surface  (all  R2C=CR2).  The  heats  of  reaction  for  R2CH-CHR2  — > 
R2C=CR2  +  H2  for  the  isolated  R2C=CR2  and  clean  surface  cases  were  also  very  similar  at  44.6 
and  46.7  kcal/mol,  respectively.  The  slightly  lower  enthalpy  of  the  isolated  R2C=CR2  may  be  due 
to  weak  interactions  between  the  double  bond  and  hydrogen  atoms  in  adjacent  unit  cells. 

Considerably  less  energy  is  required  to  remove  the  second  hydrogen  atom  on  the  dimer  than 
the  first,  as  evidenced  by  a  much  more  negative  abstraction  reaction  enthalpy,  -46.8  kcal/mol  (Eq. 
(14)).  Put  differently,  we  calculate  the  first  C-H  (R2CH-CR2-H)  bond  enthalpy  to  be  91.4  kcal/mol 
(similar  to  t-C:4H9-H,  95.8  kcal/mol^S),  while  the  second  C-H  (R2C-CR2-H)  bond  enthalpy  is 
only  57.4  kcal/mol.  For  comparison,  the  C-H  bond  enthalpies  in  the  moleculai  radicals  CH2CH2- 
H,  CH3CHCH2-H,  and  (CH3)2  C  CIH2-H  are  all  between  36.0  and  36.5  kcal/mol.  The  higher  C- 
H  bond  energy  associated  with  the  surface  species  is  mainly  due  to  the  strain  present  in  the 
abstraction  product  (R2C=CR2,  Fig.  2b)-the  surface  dimer  is  unable  to  form  as  strong  a  double 
bond  as  in  free  molecules.  Brenner  calculated  a  similar  A//  for  reaction  (14),  -50  kcal/mol  at  0  K.19 


The  implication  of  the  much  weaker  C-H  bond  in  R2CH-  C  R2  is  that  H-atom  vacancies  will 
tend  to  occur  in  pairs,  so  that  under  growth  conditions  the  (100)  surface  is  likely  to  consist 
predominantly  of  (2xl):H  monohydride  dimers  and  R2C=CR2  dimers  rather  than  isolated  radical 
sites.  One  mechanism  for  loss  of  the  second  H,  which  to  our  awareness  has  not  yet  been 
considered  in  (100)  growth  models,  is  atomic  desorption,  analogous  to  P-scission  in  molecules: 

R2CH-CR2  ->  R2C=CR2  +H  A//298  =  57.4kcal/mol  (17) 

Assuming  no  activation  barrier  for  the  reverse  reaction,  so  that  the  activation  energy  is  equal  to  the 
enthalpy  of  reaction,  and  a  preexponential  factor  of  10*^  s*^  we  estimate  a  rate  constant  of  200  s'l 
for  reaction  (17)  at  a  substrate  temperature  of  900  *C,  which  is  likely  to  be  important  Another 
likely  process  is  H  migration  to  an  adjacent  radical  site,  resulting  in: 

2R2CH-CR2  R2C=CR2  +  R2CH-CHR2  A//298  = -34.0kcal/mol  (18) 
The  driving  force  for  pairing  of  hydrogen  atoms  (and  vacancies)  on  dimers  is  a  direct  consequence 
of  the  7C-bond  in  unoccupied  dimers.  We  have  recently  noted  a  similar  tendency  for  double¬ 
occupation  of  surface  dimers  on  Si(lOO),  although  there  the  energetic  driving  force  is  much  less.28 

Turning  to  the  products  of  H  atom  addition  to  the  monohydride,  we  see  that  steric  repulsitm  of 
hydrogen  perturbs  the  structure  and  energetics  but  is  much  less  dominant  than  on  the  full  (lxl):2H 
dihydride.  The  Hr-C2-H2'  bond  angle  (101.5*)  in  R2CH2CHR2  (Fig.  2(c))  is  only  slightly  less 
than  the  unstrained  value  of  106.7*,20a  and  the  orientation  of  the  Cl-Hl  bond  is  only  slightly 
perturbed  from  that  of  an  unstrained  secondary  radical.^^  The  enthalpy  change  for  reaction  (15)  is 
slightly  more  favorable  than  that  for  the  analogous  reaction  of  2,3-dimethylbutane: 

(CH3)2CH-CH(CH3)2  +  H  (CH3)2CH2  +  (CH3)2  C  H  +  H2  A//  298  =  -1 3.4  kcal/mol  (19) 

This  implies  that  the  reduction  in  strain  energy  achieved  by  breaking  the  dimer  bond  is  slightly 
larger  than  the  increase  in  H-H  repulsion  energy. 

Steric  H-H  repulsion  is  more  important  in  the  R2CH2CH2R2  structure  (Fig.  2(d)),  causing 
considerable  distortions  in  the  C-H  bond  angles.  Also  indicative  of  steric  repulsion  in  the 
R2CH2CH2R2  surface  species  is  the  observation  that  dJi  for  reaction  (12)  is  some  20  kcal/mol 
more  positive  than  that  for 

(CH3)2CH-CH(CH3)2  +  H2  ^  2  (CH3)2CH2  AW  ^8  =  -7.15  kcal/mol  (20) 

despite  the  release  of  strain  energy  in  the  monohydride.  On  the  other  hand,  reaction  (12)  is  much 
less  unfavorable  than  reaction  (4),  forming  the  (lxl):2H  full  dihydiide,  in  which  the  H-H  steric 
repulsion  is  extreme.  Two  isolated  dihydride  units  can  reduce  their  steric  repulsion  by  simply  tilting 
away  from  one  another,  which  is  prevented  in  the  (lxl):2H  structure  by  adjacent  dihydride  units. 


Although  reactions  (15)  and  (16),  which  result  in  hydrogenolysis  of  the  dimer  bond  to  form 
two  neighboring  dihydride  species,  have  favorable  energetics,  we  suggest  that  they  may  be  minOT 
processes  under  CVD  conditions  and  that  dihydrides  will  only  exist  in  very  low  concentrations  on 
flat  (100)  terraces  (but  still  could  exist  at  steps  and  defect  sites).  Focusing  first  on  reaction  (IS), 
insertion  of  a  H  atom  into  the  dimer  bond,  a  rough  estimate  of  the  entropy  of  reaction  follows  from 
the  lost  uanslational  entropy  of  the  gas-phase  H  atom, « -27.7  caVmol  K‘^24  implying  AG*  « -6 
kcal/mol  at  298  K.  At  the  elevated  temperatures  typical  of  CVD  growth,  AG*  will  become  positive 
(cf.  Table  I),  implying  a  small  equilibrium  constant  for  reaction  (15).  An  additional  difficulty  with 
reaction  (15)  is  that  analogous  reactions  in  the  gas  phase  are  simply  not  known-abstraction 
reactions  (cf.  Eq.  (13))  dominate.  The  transition  state  associated  with  reaction  (IS)  would  be 
interesting,  most  likely  involving  a  5-coordinated  carbon  atom  or  bridge-bonded  H,  but  has  not 
been  established. 

Assuming  that  R2CH2  CHR2  can  be  produced  by  reaction  (15),  the  dihydride  could  be 
stabilized  by  a  radical  recombination  reaction  (e.g.,  Eq.  (16)).  Several  kinetic  processes  would 
compete  with  the  recombination,  however.  An  H  atom  could  migrate  to  an  adjoining  R2CH-  C  R2 
site.  Alternatively,  H2  loss  by  a-elimination  may  be  possible: 

R2CH2CHR2  R2C-CHR2  +  H2  a// 298  = -1.7  kcal/mol  (21) 

The  heat  of  reaction  for  reaction  (21)  follows  by  subtracting  reaction  (9)  from  reaction  (11).  Again, 
an  analogous  gas  phase  reaction  is  not  known,  and  the  rate  constant  cannot  be  easily  estimated. 
Taking  the  entropy  of  H2  in  the  standard  state,  29.5  cal/mol  K'1,24  as  a  rough  estimate  of  the 
reaction  entropy,  we  obtain  AG  298  *  -10.5  kcal/mol  for  reaction  (21),  indicating  a  small 
equilibrium  concentration  of  R2CH2  C  HR2. 

More-stable  dihydride  (CH2)  units  separated  by  monohydride  dimer  units  could  be  generated 
via  several  pathways.  Formation  of  R2CH2CH2R2  followed  by  a  double  migration  of  a  hydrogen 
atom  to  an  adjacent  R2CH-CHR2  dimer  seems  unlikely,  as  the  first  migration  step,  R2CH2CH2R2 
+  R2CH-CHR2  — >  2  R2CH2  C  HR2,  is  quite  unfavorable  (A//  298  =  62.6  kcal/mol).  A  more  likely 
pathway  involves  reaction  (13)  occurring  on  adjacent  dimers: 

HHHHH  HHHHHH 

Y.  ^  AH298- -98kcaVn»l  (22) 

The  estimate  of  the  heat  of  reaction  for  reaction  (22)  follows  from  the  heats  of  formation  of 
R2CH2CHR2  (Eq.  (9))  and  the  (3xl):1.33H  dihydride  (Eq.  (5))  from  the  monohydride,  and  is 
virtually  identical  to  the  negative  of  a  secondary  C-H  bond  enthalpy  .26  Formation  of  dihydride 


units  separated  by  more  than  one  monohydride  dimer  could  result  from  “migration”  of  the  radical 
center  by  rearrangement  of  C-C  dimer  bonds: 


H  H  H 


H  H 

H 

H  H  H  H  h' 

t  C- 

h—h—c 

h.  »■ 

j-K  J\  n 

H  H  H  H 
\  /  \  / 

— .  C  C-C 


H 

k 


A//298  »  Okcal/mol  (23) 


The  concentration  of  (3xl)-type  dihydride  units  will  be  determined  by  the  concentration  of 
R2CH2  C  HR2  and  by  competition  between  reaction  (22)  and  abstraction  of  H  in  the  dihydride  by 
incident  H  atoms,  but  our  estimates  of  AG*  (Table  1(c))  suggest  that  under  growth  conditions  their 
concentration  will  be  low. 


We  suggest  that  dihydride  formation  may  result  from  reactions  (15),  (22),  and  (23),  but  the 
kinetics  of  these  reactions  are  highly  uncertain  at  present.  It  is  noteworthy  that  Thomas  etcd}^  had 
difficulty  finding  evidence  for  formation  of  the  dihydride.  Further  experiments  are  clearly  required, 
and  quantum  chemical  estimates  of  the  activation  barriers  to  reactions  (IS),  (22),  and  (23)  would  be 
very  illuminating.  Finally,  accurate  frrst-principles  quantum  mechanical  calculations  of  the  structure 
and  energetics  of  selected  surface  species  are  needed  to  establish  the  degree  of  accuracy  of  MM3  as 
applied  to  surface  chemistry  on  diamond. 
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TABLES 

TABLE  I.  Thermochemical  properties  of  diamond(100)-hydrogen  reactions. 


(a)  C(100)-(2xl)  +  H2  ^ 

C(100)-(2xl);H 

T(K) 

AH’  (kcal/mol) 

AS’  (cal/mol  K-1) 

AG’  (kcal/mol) 

298 

-46.7 

-28.5 

-  38.2 

400 

-47.2 

-29.9 

-  35.2 

500 

-47.4 

-30.4 

-  32.2 

600 

-47.6 

-30.4 

-28.3 

800 

-47.5 

-30.6 

-23.1 

1000 

-47.1 

-28.6 

-  18.5 

1500 

-46.3 

-27.9 

-4.4 

(b)  C(100)-(2x1):H+  H2 

2  C(100)-(lxl):2H 

r(K) 

AH’  (kcal/mol) 

AS’  (cal/mol  K*l) 

AG’  (kcal/nx)l) 

298 

49.2 

-  11.6 

52.7 

400 

49.9 

-9.6 

53.7 

500 

50.4 

-  8.5 

54.7 

600 

50.7 

-7.9 

55.4 

800 

51.2 

-7.2 

57.0 

(c)  3  C(100)-(2xl):H  +  H2  - 

>  2  C(100)-(3xl):1.33H 

T(K) 

AH’  (kcal/mol) 

AS’  (cal/mol  K*!) 

AG’  (kcal/mol) 

298 

-  15.6 

-42.8 

-2.8 

400 

-  14.9 

-40.8 

1.4 

500 

-  14.4 

-  39.7 

5.5 

600 

-  14.1 

-39.1 

9.4 

800 

-  13.6 

-  38.4 

17.1 
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TABLE  11.  Atomic  positions  of  surface  atoms  in  modified  (l(X))-(2xl):H  monohydride. 

Displacements  of  carbon  atoms  are  given  with  respect  their  ideal  bulk-terminated  positions, 
while  hydrogen  atom  positions  are  listed  with  respect  to  the  carbon  atom  to  which  they  are 
bonded.  Atom  labels  refer  to  Fig.  2;  an  asterisk  refers  to  a  radical  center.  The  structures 
shown  in  Figs,  lb,  2b  and  2d  are  symmetric  about  the  dimer  center. 

Displacements  (A)  along  [Oil]  [Oil]  [100] 


Cl 

0.447 

(2xl):H  (Fig.  lb) 

0.000 

-  0.057 

C3 

0.066 

0.000 

0.027 

HI 

-  0.437 

0.000 

1.021 

Cl* 

0.406 

(2xl):H  -  H  (Fig.  2a) 

0.000 

-  0.020 

C2 

-  0.380 

0.000 

-  0.048 

H 

-  0.478 

0.000 

1.014 

Cl 

0.519 

(2xl):H  -  2H  (Fig.  2b) 

0.000 

-  0.106 

Cl* 

0.033 

(2xl):H  +  H  (Fig.  2c) 

0.000 

-  0.069 

C2 

0.044 

0.000 

0.007 

HI 

-  0.380 

0.000 

1.034 

HI' 

-  0.933 

0.000 

0.566 

H2' 

0.746 

0.000 

0.807 

Cl 

-  0.023 

(2xl):H  +  2H  (Fig.  2d) 

0.000 

-  0.014 

HI 

-  1.005 

0.000 

0.399 

H2 

0.540 

0.000 

0.949 

Top  view: 


[Oil] 


Side  view:  llOOJ 


L 


[Oil] 


\ci  cij 

(b)  HI  H2 


1133*  W 
1.63  1  11 


^\26.4 


\  Cl  C5l 

(d)  HI  H2  H3i|M%fll  H4 


109.4'y 

139 /C7  1  11 


FIG.  1.  Top  and  side  views  of  atomic  structures  of  clean  and  hydrogenated  diamond  (100) 
surfaces:  (a)  (2x1);  (b)  (2xl):H;  (c)  (lxl):2H;  and  (d)  (3xl):1.33H,  Bond  lengths:  A. 

O,  and  0:  carbon  atoms  in  top,  second,  and  third  layers,  respectively.  0:  hydrogen  atoms. 


